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ABSTRACT 

Recent numerical simulations indicate the presence of dynamical instabilities of the f-mode in dif- 
ferentially rotating stars even at very low values of T/|VF|, the ratio of kinetic to potential energy. In 
this Letter we argue that these may be shear instabilities which occur when the degree of differential 
rotation exceeds a critical value and the f-mode develops a corotation point associated with the pres- 
ence of a continuous spectrum. Our explanation, which is supported by detailed studies of a simple 
shell model, offers a straightforward way of understanding all of the key features of these instabilities. 
Subject headings: hydrodynamics — instabilities — gravitational waves — stars: neutron — stars: rotation 



1. INTRODUCTION 

Interest in dynamical instabilities of rotating stars is 
strong because of the prospects for detecting gravita- 
tional waves from the oscillations of neutron stars. Mod- 
elling suggests that for gravitational wave emission to be 
at a detectable level we require that the oscillations be 
unstable, and hence capable of growth. 

One key factor under investigation is the effect of dif- 
ferential rotation. Differential rotation may arise in neu- 
tron stars in several circumstances. The first is at birth; 
the latest studies of rotational core collapse indicate 
that neut ron stars will be born with strong differential 
rotat ion l|Dimmelmeier. Font fc Miilled 120021: lOtt et all 
120041) . Subsequent accretion of superno va fallback mate - 
rial or material from a companion star l|Fuiimotolll99^ . 
may also drive differential rotation, at least in the sur- 
face layers. Another possibility is that oscillations may 
drive t he star into differential rotation, via non-linea r 
effects ijRezzolla et a]|l2000l iLevin fc Ushomirskvl 120011) . 
A differentially rotating massive neutron star may also 
be generated following a binary neutron star merger 
ijShibata fc Urvuil2000(l . Differential rotation will be rel- 
evant provided that it can be maintained in the presence 
of viscosity or an internal magnetic fiel d, which will a ct 
to bring the star into uniform rotation l* ShaDiro|l200Cf) . 

Differential rotation has two key effects that are rele- 
vant to the study of dynamical instabilities. Firstly, dif- 
ferentially rotating stars exhibit a continuous spectrum, 
with dynamical behaviour that is distinct from the dis- 
crete normal modes of oscillation found in uniformly ro- 
tating stars. Secondly, differential rotation can lead to 
the occurrence of dynamical shear instabilities, unstable 
oscillations that do not exist in unifo rmly rotating sys- 
tems (Balbinskll lT985bl: lLuvte^ll990a^ . 

Before proceeding we note some points on terminology. 
We state that a mode is co-rotating if its pattern speed 
Up — cr/m is positive, and counter-rotating if its pattern 
speed is negative (where the modes are assumed to be- 
have as exp(— i(o't — rmf))). If the pattern speed of the 
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mode matches the local angular velocity at a point in the 
star we state that the mode has a corotation point. The 
range of frequencies in which modes possess corotation 
points is termed the corotation hand. All modes with 
corotation points are co-rotating, but not all co-rotating 
modes possess corotation points. 

We will also make reference to the degree of differen- 
tial rotation. This refers to the difference between the 
maximum and minimum angular velocities within a dif- 
ferentially rotating star. Near uniform rotation the dif- 
ference is small, and the degree of differential rotation is 
low. The difference increases as the degree of differential 
rotation rises. 

The majority of studies of the oscillations of differen- 
tially rotating systems have focused on the I = m = 2 
f-mode (the so-called bar mode). As in a uniformly 
rotating star, the f-mode is split into co-rotating and 
counter-rotating branches (Fig. QJ. As the ratio of ki- 
netic to potential energy, f3 = T/|W^|, increases, the 
counter-rotating mode pattern speed becomes less neg- 
ative, passing through zero at the point (3 — Ps. It is 
at this point that secular instability would arise in the 
presence of gravitational radiation reaction. In the ab- 
sence of such mechanisms, the counter-rotating branch 
pattern speed continues to increase until it merges with 
the CO- rotating branch at j3 — (3d- This merger gives 
rise to the well-studied dynamical bar mode instability 
tg^hibata. Baumgartc & S hapiro 2000 ; New & Shapi^ 
ifiOlHYoshida el 31200 jUKaririo Erigu chi 2003). The 
value of Pd is high (> 0.2) even for high differential ro- 
tation. 

Until recently, the high (3 bar mode instability 
was the main focus of attention . Re cent studies by 
iShibata. Karino fc Eriguchil l|2002L 12003(1 have identified 
dynamical bar-mode instabilities in highly differentially 
rotating Newtonian polytropes at values of (3 as low as 
~ 0.01. The mechanism behind these instabilities has 
yet to be identified. 

The presence of corotation points may be critical. In 
differentially rotating stars the co-rotating branch of the 
f-mode traverses the corotation band as [3 is increased 
(Fig. HI . In this Letter we demonstrate that the reported 
features of the low (3 instabilities can be explained if the 
co-rotating f-mode develops a dynamical shear instabil- 
ity at the point where it enters the corotation band. This 
mechanism has been studied in detail in a simple shell 
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Fig. 1. — Typical dependence of pattern speed ap on /3 for the 
I = m = 2 f-modes of a differentially rotating stellar model. Tlie 
quantity Qc is discussed in the text. Plot (a) depicts low or in- 
termediate differential rotation, plot (b) high differential rotation. 
The corotation band, traversed by the co-rotating f-mode, gets 
wider as the degree of differential rotation rises. 



mode l llWatts et a]|l2003t IWatts. Andersson fc Williamsl 
|2004() . Dynamical instabilities associated with the 
presence of corotation points are also w ell known 
in accretion disc theory LPapa l oiz^u fc Prin gle 1985; 
iGoldreich. Goodman fc NaravanI Il986|) . This would 
however be the first occasion on which they have been 
observed in stellar models. 

In section|21we explain the natm'c of oscillations within 
the corotation band, and review the results of a simple 
model in which the corotation instability mechanism op- 
erates. In section |21 we show how the mechanism can 
account for the unusual features of the low /3 bar-mode 
instabilities, using the supporting evidence provided by 
Isiiibata . Karino & Eriguchi ( 2002, 2003). 

2. THE NATURE OF OSCILLATfONS IN THE 
COROTATION BAND 



In two previous papers llWatts et all I200I 
IWatts. And ersson & Willian^ l2004() we examined 
the nature of oscillations within the corotation band 
for a differentially rotating incompressible spherical 
shell, by solving the linearized Newtonian perturbation 
equations. We found that within the band there was 
a continuous spectrum that is physically distinct from 
the more familiar discrete normal modes. In addition, 
we found dynamical shear instabilities that set in when 
modes developed corotation points. In this section 
we will show that the nature of solutions within the 
corotation band remains the same for a more realistic 
stellar model. 

We assume that the star is an inviscid polytrope 
and that the perturbations are adiabatic and non- 
axisymmetric with (^-dependence exp{imip). We also 
assume that the perturbations have time dependence 
exp(— zcrt), (T being the frequency. We define the pres- 
sure P, density p, sound speed c, gravitational potential 
$, and angular velocity f2. We work in the inertial frame 
and use cylindrical coordinates {w,ip,z). Commas are 
used to denote partial derivatives. Combining the Euler 
equations, continuity equation and a poly tropic equation 
of state for the perturbations, the Eulcrian perturba- 
tion 6(f> = SP/ p obeys the following equation ijBalbinskil 
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where the corotation parameter a = mfJ — cr, the 
Lindblad parameter L = 2^Cl — (i^, and the vorticity 
Vl = 20. + w(Kl/dw. Note that for a barotropic star 
dil/dz = dVL/dif = 0, so that the fluid rotates on cylin- 
ders Vt = Vl{vj). 

Solution of equation for real a is problematic at 
the corotation radius vj = zuc where (7 = 0, and at the 
Lindblad radii, tu = tul where L = 0. In the neigh- 
bourhood of these points the character of the solution 
can be investigated using Frobenius analysis. The solu- 
tion is regular at the Lindblad radii. At the corotation 
radius, although (j) is continuous its first radial deriva- 
tive is not. In general, there is both a finite step and 
a logarithmic discontinuity in the first derivative at Wc 
(sec Watts ct al (2003) for more detail). The additional 
degree of freedom introduced by the finite step allows us 
to construct a continuous spectrum of "eigenfunctions" 
for all frequencies within the corotation band. 

The continuous spectrum eigenfunctions found in this 
more realistic stellar model have the same mathemati- 
cal character as those found in the shell model. As in 
this simple model, if one solves the initial value problem 
one finds that the collective perturbation associated with 
the continuous spectrum is non-singular and hence phys- 
ical ijWatts. An dersson & Williams (2004), see also .Cas^ 
1|1 96011 : IBalbiTTsk i (1984)). 

The dynamical shear instabilities found in the simple 
shell model occurred where modes crossed into the coro- 
tation band, at the point where modes merged with the 
continuous spectrum. Given that the nature of the con- 
tinuous spectrum remains the same in the more realistic 
model it is legitimate to ask whether the same instabil- 
ity mechanism might also operate. In section |31 we will 
argue that the low (3 instabilities exhibit all of the key 
features of this mechanism. In order to do this we must 
first review some pertinent results from the shell model 
^Watts et a,1ll200a iWa,tts. Andersson Wi1lia,mFll2004ll. 

As the degree of differential rotation increases, many 
of the normal mode frequencies approach the edge of the 
expanding corotation band. Some approach the bound- 
ary tangentially and do not cross into corotation. Those 
that cross into corotation at low degrees of differential 
rotation remain stable. Their eigenfunctions generally 
acquire a degree of singularity, but can be distinguished 
from the rest of the continuous spectrum by the fact that 
they do not have a finite step in the first derivative. This 
led us to term these solutions "zero-step solutions" . We 
then carried out numerical time evolutions and an an- 
alytical evaluation of the initial value problem. When 
we produced power spectra of the time evolution data, 
the zero-step frequencies stood out just as clearly as the 
discrete normal modes outside the corotation band. The 
solutions that one finds when one solves the initial value 
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Fig. 2. — A schematic representation of Fig. 3 of 

IShibata. Karino &: Erig uchi ( 2003) showing the regions of parame- 
ter space occupied by the high fi bar mode and low /3 instabiUties. 
The Unes marked A, B and C are discussed in the text. 



problem are in fact non-singular, proving that these so- 
lutions have physical relevance. 

At high differential rotation zero-step solutions are still 
generated where modes cross into the corotation band. 
In addition, dynamically unstable (complex frequency) 
modes can be generated at the crossing point for certain 
rotation profiles. The eigenfunctions of these dynami- 
cally unstable modes are perfectly regular (recall that our 
earlier analysis applied only to real frequencies). Their 
pattern speeds are however in the corotation band. We 
refer to this instability mechanism as the corotation in- 
stability mechanism. 

3. DYNAMICAL INSTABILITIES 

In reviewing the low f3 instabilities we will focus on 
the j-constant rotation law used by Shibata et al. The 
j-constant rotation law is 



(2) 



The parameter A'^ sets the degree of differential rotation, 
uniform rotation being reached in the limit as >1 — > oo. 
As ^ the degree of differential rotation increases. 
For a given equilibrium model, ilc is a function of the 
parameter (3. 

Following the convention used by Shibata et al and 
working in dimensionless units where < -nj < 1, the 
corotation band occupies the range of frequencies 
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Note that the corotation band widens as de- 
gree of differential rotat io n ris es. Fig. 3 of 
IShibata. Karino Sz Eriguchil l)2003j) .which we repro- 
duce in schematic form in Fig. |21 illustrates the regions 
of parameter space occupied by the high /3 and low (3 
instabilities. The following features require explanation. 

For low degrees of differential rotation, the low (3 in- 
stabilities are not observed. Dynamical instability does 
not set in until high where the two branches of the f- 
mode merge (Track A, Fig. 12)1 . For intermediate degrees 
of differential rotation, the co-rotating f-mode becomes 
unstable at low (3. It stabilises again at intermediate (3 
before merging with the other branch of the f-mode and 



becoming dynamically unstable again at high (3 (Track 
B, Fig. 121). For the highest degrees of differential rota- 
tion there is no stable region after the onset of the low 
P instability (Track C, Fig. With this in mind, we 
propose the following explanation of these three feature. 

The fact that no instabilities are observed for stars 
with low degrees of differential rotation (Track A, Fig.[2Jl 
suggests that, as for the shell model, this type of shear 
instability can only develop if the f-mode crosses into the 
corotation band when the degree of differential rotation 
exceeds a certain threshold. For the shell model we were 
able to determine this threshold analytically. Develop- 
ing a similar criterion for the more realistic model will 
however be more difhcult. One may, in fact, have to re- 
sort to numerical calculations. Below the threshold the 
f-mode passes into the corotation band but remains sta- 
ble. Whether the eigenfunction is then purely regular, or 
possesses the singular character of the zero-step solutions 
that we found on the shell, is a matter for further study. 

At intermediate degrees of differential rotation (Track 
B, Fig. I2J) the co-rotating f-mode enters the corotation 
band and goes dynamically unstable. The lower bound 
of the low P instability region marks the point at which 
the co-rotating f-mode crosses the upper edge of the coro- 
tation band. The upper edge of the low (3 unstable region 
marks the point at which the co-rotating f-mode exits the 
corotation band. The system is then stable until the co- 
rotating branch merges with the counter-rotating branch 
at high /3 (see Fig.^a), which illustrates the passage of 
the co-rotating f-mode through the corotation band and 
its exit prior to mode merger). 

As the degree of differential rotation rises (Track C, 
Fig. 12), the corotation band gets wider. For the scaling 
used in Fig. ^ for the j-constant law, the upper edge of 
the corotation band remains fixed whilst the lower edge 
approaches the zero-axis. The counter-rotating branch 
of the f-mode will now enter the corotation band before 
the co-rotating branch can emerge (see Fig.^b)). There 
is thus no region of stability. This raises two interesting 
questions. Does the counter-rotating f-mode also develop 
dynamical instability when it enters the corotation band? 
And what becomes of the traditional bar mode instability 
at high /3? If the co-rotating f-mode never stabilises then 
the traditional mode merger route is excluded, but we 
cannot rule out the possibility of additional instabilities 
arising via zero-step mergers (Balbinski 1985b). 

The corotation mechanism provides a simple qual- 
itative explanation of the location the unstable re- 
gions of parameter space. It also leads to quantita- 
tive and qualitative predictions that we can test against 
the results of the numerical sin i ulatio ns reported by 
'Shib ata. Karino & Eriguchil ll2?ffll l200l . 

Firstly, all of the observed low /3 instabilities should 
have real parts of frequency that lie within the corota- 
tion band. Further communication with the authors of 
Shibata. Karino & Eriguchi (2002, 2003) has confirmed 
that all of the low (3 unstable modes found for the j- 
constant law lie within the corotation band (S. Karino, 
private communication) . 

To reproduce the shape of the instability region the 
co-rotating f-mode must enter the corotation band at 
lower values f3 and exit the band at higher values of (3 as 
the degree of differential rotation rises. Using numerical 
modelling we have confirmed that this is the case. 
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The third prediction is that the boundaries of the low /3 
region of parameter space represent the points at which 
the f-mode enters or leaves the corotation band. For 
models on the lower edge of the region, we should find 
that (Tr/^c = ■mA'^/{A'^ + 1), while for models on the 
upper edge of the low /3 region we should find ar/ftc = w, 
where (7^ is the real part of frequency. This prediction, 
which would provide the most clearcut confirmation of 
the corotation mechanism, is however the hardest to test. 

The problem lies in the growth times. Although the 
corotation instabilities are genuine dynamical instabili- 
ties (in that they arise in systems for which no dissipa- 
tive mechanism is present), they only exhibit dynamical 
(short) growth times deep inside the corotation band. 
The shell model was sufficiently simple that we were 
able to track unstable modes with \ai/ar\ 10~^ (at 
being the imaginary part of frequency). It was the abil- 
ity to detect these modes with very long growth times 
that allowed us to pinpoint the onset of instability as 
being near coincident with the edges of the band. Track- 
ing slow growing instabilities using the types of numer- 
ical codes used to investigate realistic stellar models is 
more difficult. The code used to generate the results of 
Shibata et al, for example, was sensitive to modes with 
Wi/o'r\ ^ 10^^ and larger (S.Karino, private communi- 
cation). Such a limitation on the shell study would have 
made it impossible for us to detect the unstable modes 
until they were well inside the band. Conclusive tests 
of the third prediction will therefore require improved 
sensitivity to modes with long growth times. Nonethe- 
less those low f3 instabilities for which growth times have 
been published do exhibit the expected behavior in that 
the growth times are shortest in the center of the band, 
ge tting longer towards the edges (see for example Fig. 4 
of lShibata. Karino fc Eriguchil l|200^ 'l. 

Although there are other inst ability mech anis ms asso- 
ciated with the corotation band (Balbinski 1985b; LuvtenI 
Il990al IPapaloizou fc Pringlail985^1 . all have hmitations 
as an explanation for the low /3 bar mode instabili- 
ties. In these scenarios the co-rotating f-mode would 
enter the corotation band as cither a zero-step or de- 
caying mode (although no decay is observed), merg- 
ing inside the band with another zero-step or decaying 
mode to go unstable. What it could be merging with 
is however unclear: the counter-rotating f-mode is out- 
side the band at low 0, and the r-modes never enter 
ijKarino. Yoshida fc Eriguchill200lT) . In addition it is not 



Balbinski E., 1984, MNRAS, 209, 145 
Balbinski E., 1985, A&A, 149, 487 
Balbinski E., 1985, MNRAS, 216, 897 
Case, K.M., 1960, Physics of Fluids, 3, 143 

Centrella J.M., New K.C.B., Lowe L.L., Brown D.J., 2001, ApJ, 
550, L193 

Dimmelmeier H., Font J. A., Miiller E., 2002,A&A, 388, 917 
Fujimoto M.Y., 1993, ApJ, 419, 768 

Goldreich P., Goodman J., Narayan R., 1987, MNRAS, 221, 339 
Karino S., Eriguchi Y., 2003, ApJ, 592, 1119 

Karino S., Yoshida S., Eriguchi Y., 2001, Phys. Rev. D, 64, 024003 

Levin Y., Ushomirsky G., 2001, MNRAS, 322, 515 

Luyten P.J., 1990, MNRAS, 242, 447 

New K.C.B., Shapiro S.L., 2001, ApJ, 548, 439 

Ott CD., Burrows A., Livne E., Walder R., 2004, ApJ, 600, 834 



clear whether the unstable mode could retain a bar-like 
character if it arose through merger with another mode 
type; the corotation mechanism, by contrast, permits 
bar-like instability as the f-mode would not merge with 
any other mode. 

The corotation mechanism may also be responsi- 
ble for the low /? m 1 spiral instabilities ob- 
se rved by several authors, since the unstable modes 
of iPickett. Duriscn 'fc DavisI l|199fiD and ICcntrcU a et al 
()200l(l are all within the corotation band. This is a topic 
for future work. 

4. CONCLUSIONS 

The presence of shear instabilities in differentially ro- 
tating stars is conceptually very interesting. It is also 
possible that these instabilities will turn out to be of 
considerable astrophysical significance. These instabili- 
ties may limit the range of astrophysical rotation laws. 
In addition the dynamically unstable modes may lead to 
interesting gravitational wave signals. It is worth noting 
that, while it is difficult to envisage scenarios that lead 
to the large value of /3 required for the more familiar bar- 
mode instability, one can easily imagine situations where 
the low P shear instabilities will operate. 

We have suggested a scenario, motivated by our anal- 
ysis of the differentially rotating shell, that explains the 
low /3 f-mode instabilities observed by Shibata and col- 
laborators. In this scenario dynamical shear instabilities 
arise when the co-rotating f-mode enters the corotation 
band and the degree of differential rotation exceeds a 
certain threshold value. The corotation mechanism pro- 
vides a straightforward explanation for all of the unusual 
features of the instabilities observed in numerical simu- 
lations, and is supported by all available data. Further 
tests of the mechanism will require high precision stud- 
ies on the onset of the instability, a task that will require 
codes sensitive to dynamically unstable modes with very 
long growth times. 
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